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NONENaATURE 
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Bohr's radius, 5.29 x 10"^ un 

total number of displacements per proton of energy E 
number of displacements per unit volume, pm'^ 
proton energy, KeV 
displacement damage threshold, eV 
Rydberg's constant, 13.6 eV 

efficiency for diffusion to junction along direction defined 
by M 

efficiency function for diffusion from depth x to the 
junction 

displacement mean free path for proton of energy E, i«i 

minority carrier diffusion length, pm 

average mass of gallium and arsenic nucleus, amu. 

nuclear density in GaAs crystal, 

distance traveled along proton path, pm 

path length of proton of energy E, pm 

minority carrier rms distance traveled before recombination, 
pm 

average projected range of a proton of energy E, pm 

energy transfered to struck nucleus, keV 

maximum energy transferred to struck nucleus, keV 

penetration depth into solar cell, pm 

average charge number of GaAs nucleus 

change in collection efficiency 

photon absorption pathlength, pm 

cosine of direction of diffusion in reaching junction 

average direction cosine of scattered proton of energy E 
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average number of displacements formed by one proton 
scattering event 

number of displace^nts per unit pathlength for a proton 
of energy E, (Mm)** 

cell collection efficiency for a normal unirradiated cell, 

displacement cross section for protons, 
recombination cross section, um^ 


A Simple Model of Proton Damage in GaAs Solar Cells 


John Vi. Wilson, Gilbert H. Walker, and R. A. Outlaw 
Langley Research Center 


SUMMARY 

A simple proton drnnage model for GaAs solar cells is derived and 
compared to experit^ntal values of change in short circuit currents. The 
recombination cross section associated with the defects was determined from 
the experimental comparison to be - 1.2 x 10"*^ cm^ in fair agreement with 
values determined from the deep level transient spectroscopy technique. 


Introduction 


Gallium arsenide solar cells have received considerable attention due 
to their potential usefulness In high-power space-energy systems as well as 
special space probe applications where high operating temperature Is a 
factor (ref. 1). However, space radiation damage to the cell may be a 
limiting factor. Consequently, radiation damage experimental studies have 
been cone ted and the present report Is an attempt to develop a simple 
model In wnich these experimental results can be understood. 

Earlier models for electron radiation damage assumed the defects to be 
produced uniformly throughout the cell volume, and modeled the cell perfor- 
mance In tsrms of cell-averaged diffusion lengths of the minority carriers 
(refs. 2 and 3). In view of the short pathlengths of the protons used In 
radiation damage studies and the specific dependence of the cell efficiency 
on the proton energy, the assumption of uniform damage within the cell 
seems unwarranted for low-energy proton exposure. Consequently, the pres- 
ent report treats geometric distribution of the displacement damage In 
detail and cell performance Is not treated by a cell -averaged parameter 
such as diffusion length. Cell performance Is evaluated In terms of the 
cell-averaged minority carrier recombination probability In diffusion to 
the junction of the cell. The average of the minority recombination proba- 
bility over the cell active region weighted according to the solar- averaged 
photoabsorption rate is used to estimate the decrement in the short circuit 
current . 
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Theory of Defect Formation 

Atomic displacements caused by proton impact with atomic nuclei result 
in crystal defects. Ute formation rate of these defects is related to 
Rutherford's cross section (ref. 4) 

a. p2 

. ,C> - ’ *0 ^2,1 1 . 
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L u m 


where do is Bohr's radius. Ep is Rydberg's constant, 1% is the 
atomic nun^er of the struck nucleus, M 2 the corresponding nuclear mass 
niM>er, E is the proton kinetic energy, Ep is the energy required to 
displace the nucleus from It's lattice site, and T,g is the maximum 
energy transfer in the collision 


4M- 

(1 + M^ 


( 2 ) 


The threshold for displacement requires T,n > Eg. The value of 
Ep - 25 eV insures that only close collisions result in displacement so 
that screening corrections to the Rutherford formula are unimportant 
(ref. 4). If the atomic recoil energy is sufficiently large (T»Ep) then 
additional displacements can be produced by the recoiling nucleus before 
coming to rest at an interstitial site. The average nutM>er of recoils 
produced by one initiating proton collision event is given as a function of 
the maximum energy transfer by 


1 + 
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log (VEp) ; 


Tm > 2^0 
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(3) 
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assuming half of the recoil energy produces further displacements and the 
other half is dissipated in other processes (ref. 4). These quantities 
allow the calculation of the number of displacements produced per unit 
distance traveled by a proton of fixed energy. 

In passing through a crystal, most of the energy of a proton is 
transferred to orbital electrons (ref. 5). The pathlength traveled in 
coming to rest is found by fitting the data of reference 5 as 

P(E) = 0.077 * 1.125 x 10"^ E^*®^ (4) 

where E is in keV and P is in urn. As derived from slowing down theory, 
a unique value of kinetic energy can be associated with each position along 
a proton's trajectory. The proton energy as a function of the distance P 
yet to be traveled before coming to rest is given by 
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E - 209.6 P^*°®/(1 ♦ 1.055 


( 5 ) 

as determined from the stopping power data of Andersen and Ziegler 
(ref. 5). In coming to rest, the proton undergoes multiple scatterings 
from atomic nuclei (of which a few result In displacements) which alters 
ever so slightly Its direction of motion. The depth of penetration R and 
pathlength P are approximately related by (ref. 5) 

P’^ R = [1 - exp (- 0.084 (6) 

This ratio Is related to the average deviation In the direction of motion 
and Is most Important at low energies. Ute average depth of penetrat1(W 
and Initial energy as related through equations (4) and (6) can be 
approximated by 

R(E) - 0.0062 E + 2.92 x 10’® E^'^^ (7) 

There Is no unique energy associated with a given depth of penetration due 
to multiple scattering. However, the average energy of protons which 
penetrate and stop at a depth R Is 

E * 593 R^*®/(1 + 3.71 R®*®) (8) 

The above quantities were used to determine the displacement density within 
a GaAs crystal. 

A proton of energy Eg incident on the face of the crystal will 
travel a distance 

Po ■ '’(y (®) 

before coming to rest. After traveling a distance p the energy will be 
reduced to 

E - 209.6 (pQ - p)^*°®/[l + 1.055 (Pq - p)^*^^] (10) 

At this position p the displacement mean free path Is 

tp(E) - 1/n 0 j,(E) (11) 

where n Is the density of scattering centers In the crystal (4.42 x 
10^71411^} and o{} Is the displacement cross section averaged for GaAs 
(M 2 « 72.5 and Z 2 * 32). The average number of displacements per unit 
pathlength Is 


Cj(E) . -V (E)/yE) 


( 12 ) 
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The use of equations (10) and (12) allow appropriate partitioning of the 
proton's energy Into electronic excitation and displacements everywhere 
along Its path. 

The number of displacements along the protons path Is rH>w to be 
related to the displacement damage In the crystal. For normal Incident 
protons of energy Eq on the face of a crystal, the number of 

displacements along Its path are^ given by equations (10) and (12). 
However, by the time its energy Is reduced to E, It has penetrated to an 
average depth x given by 

X » R(E^j) - R(E) (13) 

The path length and penetration depth are related to the average direction 
cosine (ref. 6) and Is approximated here by finding the ratio tT(E) ■ 
dP(E)/dR(E) using equations (4) and (6). In terms of m(E), the average 
number of displacements per unit depth Is 

^ . u(£) E„(E) (14) 

Where x Is found from relation (13). The effects of multiple scattering 
are demonstrated In fig. 1. The results of equation (14) for the average 
proton path 1$ shown as the full line In comparison to calculations 
neglecting multiple scattering according to equation (12) shown as the 
dashed line. 

The total number of displacements formed along the path of a proton 
with Initial energy Eq Is 

O(E^) - «n(E) ^ dE (15) 


The numerical evaluation of equation (15) Is found to be 


D(E^) 


0 

1.6 ♦ 14.5(1 ♦ 0.005 eJ'*®) logjj E^ 


£,< .728 

(16) 

Efl > .728 


where the threshold displacement energy for the proton Is 0.728 NeV. 
Equation (15) was also evaluated using the displacement theory of Lindhard, 
Scharff and Schlott (LSS) as discussed In reference 7. Tlie numerical 
values are shown In comparison to equation (16) In figure 2. 

It Is assumed that these displacements form recombination centers for 
electron>hole pairs produced by photon absorption. A pair once formed 
undergoes thermal diffusion until It recoi^lnes or Is separated at the 
junction. The relation between the rms distance traveled In moving to a 
position a distance L away from the source point Is (ref. 8) 


S 


ORiGIMAL PAGE IS 

OF POOR QUALITY 


(17) 

If Or is the recombination cross SKtIon and L the distance along an 
arbitrary straight line path to the junction, then the fractional loss of 
pairs due to recombination it, reaching the junction along a fixed direction 
Is 


f(u) 


[1 - exp(-i* 0 D (x)^ dx/ii)] 

[I - exp(. 0 ^ 0^(x)ij6 dx/u)l 


X > X, 


X < Xi 


08) 


uhere w Is the cosine of the direction to the junction. Averaging the 
fractional loss over all directions towards the junction 


F(*) • / ^(m) dll 
0 

and noting 

/‘l 0,(x)dx • ♦[0(E^) - 0(Ej)] 


results In 
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and f Is 

the 

proton 

fluenee 

D(Ej))J : 







(21) 

0(EJ)] : 

^3 




Assuming the trajectory of an average proton as defined by equation (8), 
Ex and Ej become 


E 


X 


593 (R^ - x)^^^ 

1 ♦ 3.7l(R^- x)^^^ 


( 22 ) 


and 


j 


593 (»q - 

1 ♦ 3.71 - x.)l^* 


(23) 
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CoKplete fiegitct of inuUlplt tcaUtrlng tfftcts 1$ accoMplltlied by 
replacing equations (22) and (23) by 

2.00 , 1.43, 

E - 209. 6(P - X) /[I * 1.0SS(F - x) ] (24) 

X 0 0 

and 

Ej • 209.6<Pj - + l.0SS(P5 - (25) 

If p(x) Is the nomal efficiency function for producing and collecting 
electron-hole pairs at the junction, then the efficiency change for the 
proton damaged crystal Is 

t t 

An • / p(x) F(x) dx// p(x) d(x) (26) 

0 0 

Mhere t Is the maximum depth of the active region of the cell. To 
simplify the present calculation without serious error. It Is assumed that 
p(x) Is proportional to the photoabsorption density normalized to unity as 

p{x).^*-’‘/‘ (27) 

where the solar absorption length x Is 0.714 im (ref. 9). 
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RESULTS 


The geometry of the soUr cell Is showi In figure 3. The change In 
the cell current collection efficiency as given by equations (24) and (25) 
were evaluated numerically for solar cell parameters shown. In traversing 
the cell, the protons follow neither the trajectory of the average proton 
nor the trajectory in which multiple scattering Is neglected as shown In 
figure 1. The effects of multiple scattering were estimated by averaging 
the cell damage for the two functions shown In figure 1 In which some 
effects of deviations about the average trajectory are Included. It Is 
clear that an understanding of the low-energy experimental data requires 
detailed modeling of multiple scattering effects. The winctow thickness 
parameter which varied from cell to cell in experimental tests 
(refs. 10 and 11) were assumed to be governed by a uniform distribution In 
the present calculations. The model results averaged over the window 
thickness are coa^tared to short circuit current Masurements In Irradiated 
cells in figure 4. The best value of recoad)1nal1on cross section Is found 
to be 


0^ - 1.2 X 10'^^ cm^ 
r 

which Is in fair agreement with the estimated average cross section 
- 1.06 X 10“*^ on^ determined from deep-level transient spectroscopy 
(ref. 12). 
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Fig. 3. G»As solar cell structure used In the present wdel. 
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